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[  i  ]  F  ro  m  an  ana  l  y  si  s  o  f  geo  m  agn et i  c  and  so  I  ar  w  i nd  d  ala , 
Lockwood  et  al.  [1999]  (hereinafter  refcucd  to  as  LSW99) 
reported  that  the  solar  coronal  magnetic  field  had  increased 
by  more  than  a  factor  of  two  during  the  last  century.  If  true, 
this  would  be  an  important  discovery.  Recently,  Svalgaard 
and  cliver  [2005]  (hereinafter  referred  to  as  SC05)  reported 
an  analysis  based  on  our  newly  developed  mterdiumal 
variability  (IDV)  index  of  geomagnetic  activity  which 
indicated  that  cycle  averages  of  the  solar  field  varied  no 
more  than  ~25%  over  the  same  time  interval  and  arc  now 
decreasing.  Here,  we  answer  the  criticisms  of  Lockwood  et 
al  [2006]  (hereinafter  referred  to  as  LRFS06)  to  our  paper- 
In  sum,  wc  find  their  objections  without  merit.  If  our 
prediction  that  the  next  solar  cycle  will  be  the  smallest  in 
100  years  [Svalgaard  et  a l ,,  2005]  bears  out,  this  debate 
may  be  settled  by  direct  solar  wind  measurements  within  the 
next  few  years.  In  the  following  sections  we  respond  to  the 
various  points  raised  by  LRFS06;  percentage  change,  B, 
versus  B ,  regression  technique  (including  the  effect  of 
missing  data),  and  analysis  procedure. 

1-  Long-Term  Variation  of  IMF  B 

[2]  Our  IDV  paper  [SC05]  showed  that  B  follows  the 
sunspot  number  (or  rather  a  constant  plus  the  square  root  of 
l he  sunspot  number).  The  key  new  result  of  the  LSW99 
paper  was  the  finding  that  the  Surfs  coronal  magnetic  field 
hud  an  underlying  component  that  increased  independently 
of  tile  sunspot  number  Parker  [1999,  p.  416]  commented 
on  this  point  in  the  same  issue  in  which  LSW99  appeared: 

Lockwood,  Stamper  and  Wild  use  records  of  geomagnetic 
activity  [. , .]  to  show  that  the  weak  (-0,5  millitesla)  general 
magnetic  field  of  the  Sun  has  more  than  doubled  over  the 
past  100  years.  [...],  The  new  finding  complements  the 
well-known  fact  that  the  number  of  sunspots  [., .]  has  also 
doubled  over  the  same  period  of  time.  The  general  field 
varies  much  less  with  the  )  I -year  sunspot  cycle  than  does 
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the  number  of  sunspots,  and  appears  to  have  a  different 
origin." 

[j]  The  basic  finding  of  LS W99  would  be  of  fundamental 
importance,  if  substantiated.  It  is  illustrated  in  Figure  I. 
There  are  other  features  of  Figure  1  to  which  we  shall  return 
but  the  point  which  we  wish  to  emphasize  now  is  that  for  a 
fixed  sunspot  number  (say  40)  on  the  rising  branch  of  all 
solar  cycles  during  the  20th  century,  the  quantity 
calculated  by  LSW99  exhibited  a  steady  increase.  Figure  1 
shows  the  crux  of  our  disagreement  with  LSW99  In  a 
preliminary  report,  Svalgaard  et  ai  [2003,  p,  15]  noted  that 
there  was  no  long-term  trend  in  B  "other  than  a  general 
correlation  with  the  sunspot  number/'  The  absence  of  an 
underlying  long-term  trend  in  B  implies  the  absence  of  such 
a  trend  in  FSl  contrary  to  the  central  result  of  LSW99.  The 
sunspot  correlation  with  B  is  no  surprise  [e.g.,  Wung  and 
Sheeley,  2003]. 

[4]  In  SCO  5  we  noted  that  the  temporal  evolution  of  IMF 
B  near  the  Earth  seemed  to  exhibit  the  —100  year  Gleissbcrg 
cycle  often  seen  in  the  sunspot  number  (see  Li  et  al  [2005*j 
for  a  recent  analysts).  This  is  not  surprising  because  most 
(>70%,  see  Figure  8  of  SC05)  of  the  variation  of  B  is 
explained  by  the  variation  of  die  number  of  sunspots  and 
their  magnetic  fields.  Let  us  write  the  long-term  variation  of 
B  with  time  t  as  B  ( t )  =  M  +  A  cos  (2it  t  fQ  +  ^)  where  A  is 
the  amplitude  of  the  ^100  year  cycle,  jit  its  frequency,  ip  its 
phase,  and  M  is  the  overall  mean  on  which  the  Glcissberg 
cycle  rides.  A  conservative  estimate  of  the  amplitude  is  half 
the  range  from  the  mean  of  cycle  14  to  the  mean  of  cycle 
19  yielding  the  relative  amplitude  of  the  ^100  year  cycle  of 
AIM  =  0.l3t  for  a  total  range  of  -25%  as  quoted  in  SC05. 

M  We  contrast  this  view  with  that  of  LRFS06  who 
interpret  their  result  as  a  secular  increase,  although  the 
1 1 00  years'  of  LSW99  has  now  become  —50  years, 
perhaps  reflecting  the  fact  that  direct  measurements  of  the 
solar  magnetic  field  do  not  show  any  increase  in  the 
calculated  source  surface  flux  since  1974  [Arge  et  al ., 
2002]  nor  in  the  solar  mean  field  since  1968  [Kotov  and 
Kotova ,  200 1 J.  LRFS06  calculate  an  1 1-year  running  mean 
to  quantify  the  secular  increase,  but  a  running  mean  is  too 
simple  a  filter  to  completely  remove  the  sunspot  cycle. 
Because  of  remaining  cycle  signal,  t  (including  some  plain 
noise),  the  difference  between  the  cherry-picked  absolute 
maximum  of  the  running  mean  (in  1956)  and  the  absolute 
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figure  I.  Quantity  ft  calculated  by  LSW99  from  {cm)  and  {/)  and  inferred  from  interplanetary 
measurements  {heavy  blue  line);  adapted  from  LSW99.  The  red  dots  in  the  main  plot  highlight  years 
when  the  sunspot  number  Rz  was  near  40  on  the  ascending  branch  of  the  cycles  (away  from  recurrent 
high-speed  streams).  ft  shows  a  marked  increase  independent  of  the  sunspot  number  while  B  derived 
from  IDV  does  not,  as  shown  by  the  inset  {see  text).  Discrepancies  between  the  calculated  and  observed 
flux  arc  marked  with  the  red  oval 


minimum  (in  1903)  overestimates  the  long-term  variation. 
Dividing  this  difference  by  the  minimum  value  (as  done  in 
LRFS06)  calculates  this  curious  quantity 

(lM  +  (-4  +  e™,))  -  W  ~  (A  +  emi„)J}/[M  -  (A  +  eniinJj 
^  [2  A  A  Ed]**  T  £mm)/(Af  —  A  Sfliiu) 

that,  with  unknown  c,  docs  not  seem  to  have  a  useful, 
simple  interpretation. 

2*  Estimating  the  Radial  Component, 
of  IMF  B 

M  LRES06  fault  SC05  for  misrepresenting  the  results  of 
LSW99  by  conflating  B  with  ft,  the  radial  component  of  the 
IMF  in  the  ecliptic  plane  at  the  Earth,  which  is  assumed  (by 
LSW99)  to  be  proportional  to  the  total  coronal  source  flux. 
We  counter  that  there  is  no  confusion  on  this  point,  as 
LSW99  themselves  assume  Br  to  be  proportional  to  B. 

M  Physically,  IDV  does  not  depend  on  ft,  but  on  the 
southward  component,  BXf  in  a  coordinate  system  where 
seen  from  the  Sun  the  north-south  plane  contains  the  Earth's 
dipole  axis,  or  more  precisely  on  B  (draped  around  the 
"nose”  of  the  magnetosphere)  times  some  function  of  the 
angle  between  the  Earth's  field  and  the  draped  IMF.  By 
positing  a  constant  ratio  (called  sB  by  LSW99,  who  found 
.v,j  =  0.56  based  un  solar  wind  data)  between  ft  and  ft  you 
can  “ transfer’1  that  dependency  from  B  to  ft.  It  does 
therefore  not  make  sense  to  regress  IDV  against  ft  directly 


as  done  by  LRFS06,  However,  let  us  follow  their  lead  and 
see  where  it  takes  us. 

[s]  As  LRFS06  point  out,  magnetic  clouds,  coronal  mass 
ejections,  stream-stream  interactions,  and  the  generally 
higher  variability  associated  with  solar  activity  will  lend 
to  cancel  some  of  the  radial  flux  so  that  |ft|  for  such 
conditions  will  be  underestimated.  This  may  be  seen  in 
Figure  2,  where  we  have  plotted  (red  curve)  the  ratio,  L , 
between  (yearly  means  of)  \Br\  and  B  (both  first  averaged 
over  the  1-day  interval  advocated  by  LRFS06)  as  observed 
by  spacecraft.  There  is  the  expected  clear  solar  cycle 
dependence  in  the  ratio,  L  being  smaller  at  solar  maximum 
and  larger  at  solar  minimum,  approximated  by  a  simple 
expression  (black  curve):  L  =  0,53  -  RZI20QQ<  where  Rz  is 
the  International  Sunspot  Number  on  the  Ziirich  scale. 
Under  the  assumption  that  the  processes  that  cause  cancel¬ 
lation  of  | ft. |  operated  similarly  1 00  years  ago  we  extend  the 
black  curve  back  in  time.  Also  shown  in  Figure  2  {blue 
curve)  is  the  ratio,  L  (LRFS06)  =  |ft|/ft  calculated  ftom  the 
correlations  of  B  and  |ft.|  for  IDV  derived  by  LRFSG6  using 
their  preferred  method  B  and  given  in  their  Tabic  I .  Not 
only  does  Z,(LRFSG6)  vary  oppositely  to  L{ observed)  during 
the  spacecraft  era,  but  it  also  approaches  the  physically 
unreasonably  low  value  of  0.2  during  the  very  quiet  years  at 
the  beginning  of  the  20th  century  when  there  were  almost 
none  of  the  factors  present  that  lead  to  cancellation  of  flux. 
Wc  therefore  do  not  find  any  of  the  conclusions  based  of  the 
calculation  of  fft.|  by  LRFS06  to  be  valid.  The  IDV  index 
stands  on  its  own,  has  a  simple  definition,  and  does  not  need 
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Figui  e  2,  Ratio,  L,  between  (yearly  means  of)  |£r|  and  B  (both  first  averaged  over  a  I -day  interval)  as 
observed  by  spacecraft  (red  curve),  extrapolated  from  its  observed  solar  cycle  dependence  (black  curve), 
and  calculated  from  the  predictive  correlations  of  B  and  I Br\  derived  by  LRFSQ6  from  /DFusine  their 
method  B.  '  b 


to  be  “predicted “  (LRFSG6  Figure  11)  or  corrected  using 
proxies  for  Br  or  the  aa  index. 

3.  Proper  Regression  Method 

M  LRFS06  state  that  B  is  better  determined  than  IDV and 
therefore  should  be  used  as  the  independent  variable  “with 
no  error.”  This  does  not  take  into  account  that  33%  of  the 
hourly  IMF  values  are  missing,  for  some  years  even  in 
excess  of  twice  that.  We  can  simulate  the  effect  of  missing 
data  by  using  the  coverage  pattern  for  the  13  years  with 
least  data  (on  average  42%  covered)  (excluding  1963- 
1964,  where  there  was  very  little  data)  as  a  mask  over  the 
13  years  with  the  most  data  (on  average  98%  covered)  to 
select  data  to  throw  away  before  calculating  yearly  means. 
We  then  calculate  the  RMS  variation  of  the  differences 
between  the  observed  high -coverage  years  and  simulated 
same  years  with  coverage  degraded  according  to  the  mask. 
Varying  at  random  which  years  in  both  groups  correspond, 
i  c-*  which  low- coverage  year  to  use  as  a  mask  for  which 
high-coverage  year,  we  get  several  such  RMS  values, 
clustering  around  0.24  nT  (or  3.5%  of  B).  Wc  take  this  as 
an  estimate  of  the  uncertainty  of  B  as  representative  of 
the  whole  year  for  years  with  poor  coverage  (not,  of  course, 
of  the  inherent  measurement  or  intcrspacccraft  cross- 
calibration  error  of  B).  According  to  SCG5  the  standard 
error  of  a  yearly  mean  oUDVis  0.20  nTor  2,0%  of  IDV  It  is 
therefore  clear  that  the  errors  are  comparable,  biking  Into 
account  that  some  ol  the  time  B  is  well  determined.  This 
fact  might  argue  for  using  method  C,  except  that  wc  arc  not 
exploring  the  functional  relationship,  but  merely  seek  to 
predict  B  from  IDV 

M  Because  of  the  missing  IMF  data  we  adopted  a  very 
careful  procedure  for  calculating  the  means  and  gave  a  table 
(Table  3  in  SC05)  will]  the  results,  LRFS06  state  that  SC05 
used  a  recurrence  method  to  till  the  gaps  in  the  IMF  data 
series  and  claim  that  this  is  not  a  valid  technique  because  the 
autocorrelation  function  (for  I -day  values)  of  the  IMF  is 
only  0.20  at  a  lag  of  27  days.  In  fact,  as  clearly  stated  in 
SC05,  wc  interpolated  rotation  averages  for  the  few  rota¬ 
tions  with  no  data.  For  rotation  averages  the  autocorrelation 


coefficient  is  0.69  (for  a  lag  of  one  rotation),  high  enough  to 
justify  our  procedure  for  obtaining  yearly  averages  utilizing 
the  high  positive  conservation  of  the  IMF  data.  The  auto¬ 
correlation  function  for  l -day  values  is  not  relevant  for  this 
procedure.  We  have  to  confess  that  Table  3  in  SC05  was 
constructed  with  a  slight  amount  of  additional  data  and  that 
our  Figure  6  and  equation  (2)  were  not  updated  to  include 
the  very  latest  data.  The  regression  equation  using  the  latest 
data  up  to  the  present,  B  =  (0.371  ±  0.033)  IDV  +  (2.92  ± 
0.34),  is  still  within  the  error  bounds  of  equation  (1)  of 
SC05,  B  =  (0361  ±  0,037)  IDV  +  (3,04  ±  035).  For 
predictive  purposes  Method  A  is  the  proper  regression 
method  to  use.  The  goal  is  to  best  predict  or  estimate  the 
unknown  B  from  the  known  IDV  where  the  measure  of  the 
goodness  of  the  prediction  is  how  welt  (in  the  least  squares 
sense)  the  empirical  relation  reproduces  observed  values  of 
B >  realizing  that  the  relation  is  complex  and  involves  other 
variables  as  well 

[n]  The  discussion  of  the  am  index  showing  that  the 
shapes  of  the  distributions  for  all  days  and  for  days  with 
IMF  data  are  nearly  identical  (LRFS06  Figure  3)  does  not 
demonstrate  that  there  is  a  lack  of  bias  due  to  missing  IMF 
data  for  the  single  year  averages  with  which  we  arc 
concerned,  only  that  over  40  years  the  biases  arc  diluted 
out.  In  fact,  LRFS06  Figure  4  shows  that  for  some  years  the 
difference  between  (am)  and  (gjtj)w  (tor  intervals  wdth 
coincident  IMF  data)  can  reach  10%. 

[i2]  LRFS06  goes  into  great  detail  about  the  effect  of 
outliers.  In  SC05  we  elected  not  to  remove  data  points  that 
did  not  “fit/1  By  repeatedly  removing  data  points  that  do 
not  fit  one  can  get  the  correlation  as  good  as  desired  while 
eventually  decreasing  its  statistical  significance.  From 
Figure  3  and  from  LRFS06  (e.g.,  Figure  9)  it  is  clear  that 
removing  the  five  points  with  the  greatest  deviation  between 
htted  and  observed  values  of  B  results  in  a  markedly 
improved  correlation  (R2  goes  from  0,75  to  0.89).  We  find 
for  the  years  1965-2006  (first  7  months)  using  ordinary  least 
squares  regression  of  B  on  IDV: 

(0.445  ±0.026)/£>F  +  (2.21  ±0.26)  (R"  0H9| 
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Figure  3.  Identification  of  five  “outliers"  (marked  with  green  dots)  for  the  correlation  between  IMF  B 
(red)  and  ID V.  The  regression  fit  (blue  curve)  was  made  omitting  die  outliers  as  advocated  by  LRFS06. 
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B  =  0.445  IDV  +  2.21 
R2  =  0.89 


Hi  is  result  is  very  close  to  correlation  three  in  Table  1  of 
LRFS06,  obtained  using  the  “least  median  of  squares” 
metliod  that  is  much  less  sensitive  to  outliers  (equivalent  to 
largely  ignoring  them).  Whether  or  not  it  is  permissible  to 
remove  outliers  or  greatly  reduce  their  leverage  is  a  different 
debate.  Here  wc  shall  investigate  the  effect  of  using  the 
above  regression  (or  the  almost  identical  one  calculated  by 
LRSF06).  Figure  4  shows  the  result.  Because  the  higher 
slope  compared  with  the  SC05  correlation  is  largely 
compensated  by  the  lower  offset,  the  result  is  almost 
identical  lo  our  Figure  7  in  SC05,  We  approximate  the  long* 
tenn  variation  by  a  fourth-order  polynomial  shown  by  the 
red  curve.  The  relative  amplitude,  A/M,  determined  from 
this  curve  ts  0.1  13  corresponding  to  a  total  range  of  ^23%. 
Another  fourth-order  polynomial  (black  curve)  is  fitted  to 
the  data  points  for  the  sunspot  minima  only  and  shows  a 
similar  variation  simply  downshifted  to  match  the  lower 
values  at  minima.  So,  removing  the  outliers  or  downplaying 
their  relevance  by  using  a  more  “robust”  regression 
technique  does  not  change  the  long-term  variation  reported 
m  SC 05.  LRFS06  note  that  the  outliers  do  not  appear  in 
their  calculation  of  the  interhour  variability  {IHV)  index. 
They  are  confusing  the  IHV  index  with  the  (very  different) 
ID l  index,  IH\\  am,  and  aa  depend  on  solar  wind  speed 
squared  while  the  IDV  index  does  not  vary  with  solar  wind 
speed. 

[13]  To  gauge  the  importance  of  outliers,  we  compare 
yearly  averages  of  IMF  B  inferred  in  several  different  ways. 
Figure  5  shows  the  result.  The  red  curve  shows  the 
magnitude  of  IMF  if  observed  by  spacecraft  in  the  ecliptic 
plane  near  the  earth.  The  blue  curve  shows  B  calculated 
from  IDV  using  the  regression  derived  by  removing  the  five 
outliers  (or  the  almost  identical  LMS  correlation  of 
LRFSQ6).  Blue  open  circles  show'  B  calculated  from  IDV 
without  removing  the  outliers.  The  purple  curve  shows  B 
calculated  from  the  mean  magnetic  field  of  the  Sun  as  a  star 
[Kotov  et  aL #  2002;  Sva  [guard  et  al.y  2003).  The  green  curve 
shows  B  calculated  from  negative  only  values  of  the  Dst 
index  [see  SL05],  which  is  what  IDV  really  is  a  proxy  for. 
All  of  these  curves  show  occasional  disagreements  or 
“outliers.”  Instead  of  being  classified  as  “errors”  that 


should  be  removed  or  suppressed  using  “robust”  methods, 
outliers  tell  us  something  of  the  processes  that  are  behind 
the  correlations  and  should  be  studied  in  their  own  right.  For 
example  we  note  that  three  of  the  five  outliers  occur  for 
years  where  the  IMF  coverage  was  poor  (less  than  50%).  If 
Johannes  Kepler  had  used  robust  methods  in  dealing  with 
the  8  are  min  outliers  in  Tycho  Brahe's  Mars  observations, 
he  might  not  have  discovered  that  the  p lands  revolve  in 
ellipses  rather  than  in  circles  [Monitor  and  Takemoio ,  2005]. 

4.  Independent  Confirmation  From  Dst  index 

[14]  SCO  5  noted  and  demonstrated  that  yearly  averages  of 
the  Dst  index  calculated  using  negative  values  only  is 
strongly  correlated  (R2  =  0.89)  with  the  IDV  index  and  thus 
also  with  IMF  B  (R"  —  0.72).  J.  J.  Love  (personal  commu¬ 
nication,  2006)  has  recently  derived  Dst  back  to  1905. 


Year 

Figure  4.  IMF  B  inferred  from  ID  Fusing  the  regression 
with  outliers  omitted*  Fourth -order  polynomials  tilted  to  ail 
data  (red  curve)  and  sunspot  minima  data  (open  circles) 
only  (black  curve).  The  I  l*year  running  mean  (grey  curve) 
still  shows  a  wrcak  solar  cycle  variation. 
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rigure  5.  Several  estimates  of  IMF  B:  Observed  by  spacecraft  (red).  Calculated  from  IDV  omitting 
outliers  (blue  curve),  not  omitting  outliers  (blue  circles).  Calculated  from  the  Sun-as-a-star  mean 
magnetic  field  (purple).  Calculated  from  yearly  means  of  negative  Dst  index  values  (green).  Note  that  the 
Dst  curve  does  not  have  the  outliers  in  2000-2001. 


Figure  6  shows  the  substantial  agreement  between  B  in¬ 
ferred  from  IDV  and  from  Dst  (<0).  The  correlations 
between  IDV  and  B  and  between  Dst  (<0)  and  B  have 
different  outliers  and  yet  give  substantially  the  same  in¬ 
ferred  B  over  the  last  100  years*  showing  that  with  enough 
data  straightforward  analyses  work  well. 

5,  Residuals  Versus  Observed  Values 

[is]  Figure  7  shows  that  the  residuals,  Z?„bs  —  Bt-iu 
obtained  from  regression  equation  (2)  in  SC05,  show  no 
trend.  We  follow  the  standard  (and  correct)  practice  of 
plotting  residuals  (observed  minus  fitted)  versus  the  fitted 
values.  LRFS06  (Figures  6  and  10)  also  plot  residuals 
versus  observed  values  and  make  the  statement  “[n]ote 
how  strong  the  bias  is  in  SC05*s  plot  [although  it  is 
LRFS06’s]:  the  slope  of  the  LMS  regression  fit  in  Figure  6 
(right)  [residuals  versus  observed  values]  is  jr  =  0.679  [...]. 
If  we  use  this  regression  to  correct  SCOS’s  extrapolated 
values  [. . .],  this  yields  [. . .]  by  equation  (1)  a  X  of  192%. 
[’  1  [TJhis  does  serve  to  show  the  extreme  sensitivity  of  X 
to  uncertainties  and/or  inadequacies  in  SC05’s  regression 
procedure.”  The  inadequacy  here  is  LRFS06’s  as  it  is 


inappropriate  to  plot  residuals  against  observed  values,  as 
we  show  by  the  following  precise  argument. 

[i&]  Consider  the  statistical  model  (using  standard  nota¬ 
tion  [e . g. ,  Draper  ami  Sm ith ,  1998]): 

E  (Y)  =  cX 

Cov{V)  =  a2l 

Then,  provided  that  X*X  is  nonsingular,  the  least  squares 
estimate  of  c  is 

c  =  (XTX)  rx'Y 

The  fitted  Y  is 

y  =  \c  —  x{xTx)  'xtf  =  cy 

where  C  =  X  (X'xr1  X1  is  symmetric  idernpotent,  that  is, 
CC1  =  C.  The  residual  vector  is 

F-  Y  =  (I-C)r. 


1 0  IMF  8  interred  from  IDV  and  Pel 


B  =  2.92  +  0.371  10  V 
8  IMF  obs 


18,  u  1&8D  I89C'  1900  1910  1920  1930  1940  1950  1960  1970  1900  1990  2000 


I  iguie  6,  Magnitude  B  of  the  interplanetary  magnetic  field  near  the  earth  observed  by  spacecraft  (red 
curve)  and  inferred  from  the  IDV  index  (blue  curve  and  regression  formula).  The  green  curve  shows  B 
calculated  from  an  extension  back  to  1905  of  the  ZXsf  index  (J.  J.  Love,  personal  communication,  2006) 
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Figure  7*  Residuals.  B0 ^  of  observed  yearly  means  of  IMF  B  from  values  predicted  by  the 

regression  equation  giving  B  as  a  function  of IDV (method  A)  from  SC05  plotted  against  the  fitted  values 
of  B.  Also  shown  is  the  seven-point  running  mean.  No  "outliers”  have  been  removed. 
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The  coefficient  of  detenni nation  R:  for  e  versus  Y  (and 
similarly  for  F)  is  defined  as 

RV  Y)  [Cov(K.e)/Cov(T,  T}][Cov( Y, e)/Cov(e, <}] 
f  17]  Now, 

Cov( F, e)  =  Co v(CF,  (1  -  C)Y)  «=  a1  (C  -  CCT)  =  0 

Therefore  R2  (e,  Y)  «  0  and  F,  Le.,  the  fitted  Y  is  not 
correlated  with  the  residuals  e.  On  the  other  hand, 

Cov(Y.e)  =  Cov( F,  (I  -  C)F)  =  tr{l  -  C)  ^  0 

1  here  fore  R"  (e,  Y)  /  0  and  observed  Y  is  correlated  with 
the  residuals.  As  a  curiosum  we  note  that  the  slope  of  the 
regression  line  through  the  residuals  plotted  against 
observed  values  equals  the  coefficient  of  determination. 

[Hi]  A  non  mathematical  argument:  The  residuals  arc 
uncorrelaied  with  the  fitted  values;  otherwise,  they  would 
not  be  the  residuals  of  a  least  squares  fit.  The  observed 
values  arc  the  (fitted  plus  residual)  values;  therefore  the 
observed  and  residual  values  are  correlated  and  it  is  not 
appropriate  to  plot  residuals  against  observed  values  and 
then  to  claim  that  this  correlation  is  an  inadequacy  of  our 
analysis.  In  LRFS06'S  Figure  6  it  is  clamed  that  die  fit  is 
that  “given  by  SC05.”  They  tail  to  mention  that  they  do  not 
use  the  yearly  averages  for  B  found  in  our  Table  3  but 
instead  use  their  own  values  for  (5),  A  small  matter,  in  the 
caption  lor  their  Figure  6  the  intercepts  quoted  are  switched 
around.  We  shall  let  the  reader  judge  the  difference  in 
homoscedadty  between  Figures  6  and  10  of  LRFS06  while 
noting  that  if  the  relative  errors  of  a  fit  arc  approximately 
constant,  which  is  typical  for  a  relation  between  a  driver  and 
the  ensuing  response,  the  absolute  errors  automatically 
show  hcteroscedaetty. 


6.  Physics-Based  Versus  Empirical  Approach 

l|l>J  The  "physics-based”  relationships  used  by  LSW99 
to  calculate  their  quantity,  F3,  have  at  their  core  an  ad  hoc 


parameter,  their /p  =  sr  ((/)  (aa)  4 95 )  0,363  +  cf,  where  /  is  the 
Sargent  recurrence  index  [Sat'gent,  1986],  aa  is  the  aa 
index,  angle  brackets  denote  yearly  average  values,  and  jr 
and  are  two  empirical  constants.  The  variation  of  Jp  >s 
dominated  by  that  m  die  solar  wind  speed.  Because  of  the 
constant  cF,  the  solar  wind  speed  determined  by  LSW99  is 
bounded  below  at  334  km/s  (for  (/)  =  0),  Should  the  actual 
solar  wind  speed  fall  below  that,  the  result  would  be  an 
artificial  increase  in  B  inferred  from  their  method.  LSW99 
state  that  the  primary  justification  for  adopting  this  partic¬ 
ular  relationship  was  that  it  produces  a  good  correlation  for 
the  interval  1964  1994.  We  have  calculated  J  following 
Sargent's  prescription  using  27-day  Bands  rotations.  The 
cross  correlation  is  often  negative  (23%  of  the  time)  and  the 
yearly  mean  (J)  is  negative  in  1956,  1969,  and  1992.  It  is 
not  entirely  clear  what  is  meant  by  {/).  LSW99  (p.  437) 
state  that  they  use  "Annual  means  of  Sargent's  recurrence 
index,  {/)  (defined  for  the  y'th  27-day  Carrington  [sic] 
rotation  period  as  /,  =  (I/I 3)  from  k  =  -6  to 

k  —  +6,  where  c  is  the  correlation  coefficient  between  two 
consecutive  intervals  of  twelve-hourly  ua  values)."  Wc 
calculated  {/}  over  a  year  using  the  13  rotations  prior  to 
the  rotation  whose  first  day  belongs  to  the  following  year, 
but  many  other  interpretations  arc  possible,  o.g..  to  choose 
as  thcy'th  rotation  the  one  that  contains  the  183rd  day  of  the 
year,  or  even  calculating  the  yearly  mean  of  the  thirteen  or 
fourteen  13-rotation  running  mean  values  within  the  year. 
We  do  not  know  how  LSW99  calculate  (/}0263  for  negative 
{/),  nor  do  we  reproduce  their  (/)  values  (we  reproduce 
Sargent's  values  as  precisely  as  his  published  ligurc  allows 
us  to  determine). 

[20]  To  illustrate  the  danger  of  applying  ad  hoc  relations 
outside  of  the  domain  on  which  they  are  defined,  wc  may 
note  that  the  year  2003  had  some  of  the  most  extensive 
coronal  holes  and  high-speed  streams  on  record  ((V)2o oj  = 
545  km/s,  the  highest  yearly  mean  ever  measured),  yet  (/) 
was  only  0,055  versus  the  ^0*4(1  usually  found  for  years 
with  strong  high-speed  streams,  leading  to  an  underestimate 
(using  the  method  described  in  LSW99)  of  Fby  ^70  km/s, 
and  an  overestimate  of  B  by  3  nT  (40%),  a  much  larger  error 
than  even  our  largest  residual  ( 1  nT).  We  take  this  as  an 
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I  igu  rc  8,  Effect  of  changing  the  yearly  average  {/}  from  0.40  to  0  for  the  level  of  activity  {an)  =  36.3  of 
the  year  2003  on  B  calculated  using  the  ad  hoc  fp  function  of  LSW99.  The  arrow  shows  the  observed 
value  of  {5} 2(io3 


indication  that  the  LSW99  procedure  is  not  inherently 
superior  as  claimed  by  LRFS06. 

[21]  With  {/)  as  small  as  G.055t  it  could  easily  by  chance 
have  been  even  nearer  to  zero  {like  {/)1956  -  -0,002  and 
(Oiwk  =  0.009).  It  is  instructive  to  investigate  the  effect  of 
(/)  changing  from  0.40  to  0  for  the  level  of  activity  of  the 
year  2003  shown  in  Figure  8.  The  arrow'  shows  the  observed 
value  of  (iOatKO,  which  is  incidentally  what  we  would 
expect  for  the  {/}  =  0.4  that  is  typical  for  this  phase  of  a 
declining  sunspot  cycle.  Note  how  strongly  B  grows  as  (/} 
approaches  zero  from  0.4,  from  7,9  to  2 1 .2,  and  that  most  of 
the  change  takes  place  very  close  to  zero,  c,g.,  from  13,2  to 
21,2  as  (/)  changes  from  0.0 1  to  0.00.  Such  extreme 
sensitivity  is  not  the  behavior  of  a  real  physical  system 
but  comes  from  wrapping  what  the  authors  considered  to  be 
a  physics-based  procedure  around  the  core  of  an  empirical, 
ad  hoc  relationship.  LSW99  (p,  439)  prudently  stated  that 
"we  assume  that  all  three  correlations  were  valid  at  all 
times  "  We  now  know  that  ai  least,  fp  fails  the  test  of  time, 
although  its  nonphysical  extreme  sensitivity  for  very  small 
values  of  {/}  should  have  been  dear  from  the  outset.  Figure  3 
of  LSW99  should  have  given  a  warning  signal  as  shown  in 
figure  I  The  rather  glaring  discrepancies  (highlighted  in  the 
red  oval)  between  calculated  by  LSW99,  during  1967- 
1970  and  that  calculated  from  die  IMF  measured  by  space¬ 
craft  (heavy  blue  curve)  are  caused  by  very  small  values  of 
{/)  for  that  interval  showing  failure  of  the  ad  hoc  relation¬ 
ship,  despite  the  stated  99.999999999996%  significance 
level  ol  the  correlation.  There  is  a  similar  problem  with 
1992,  for  which  we  calculate  {/)  =  -0,12.  If  we  set  {/)  =  0 
tor  this  year  (what  else  makes  sense?),  we  find,  using  the 
method  of  LSW99*  (8)|$g2  =  13.5  nT,  which  is  64%  larger 
than  the  observed  value  of  8.26  nT. 

[22]  When  we  first  calculated  (/Jiooj  and  found  it  to  be  so 
low,  we  thought  that  there  was  a  problem  with  the  act  index 
data,  so  wc  redid  the  calculation  using  the  um  index  and 
confirmed  that  there  is  no  problem  with  the  data.  The  reader 
may  wish  to  compare  our  Figure  9  with  Figure  1  of  LSW99 
to  gauge  some  of  the  differences  between  {/)  and  the  values 
calculated  by  LSW99  or  at  least  shown  in  their  figure. 
Because  of  the  large  variation  of  the  cross  correlation  from 


rotation  to  rotation,  the  effect  of  {/)  when  (/)  is  dose  to  zero 
is  very  sensitive  to  precisely  how  the  annual  means  are 
calculated,  as  should  be  clear  from  Figure  9  that  shows  the 
cross  correlation  for  each  rotation  pair  and  their  13-rotation 
running  mean. 

7.  Conclusion 

[22]  LSW99  reported  that  the  Sun's  coronal  magnetic 
field  more  than  doubled  during  the  last  century.  Our 
calculation  in  SCQ5  of  B  from  IDV  indicates  no  such 
increase  in  B ,  and  therefore  none  in  Bt  and  the  coronal 
magnetic  field.  The  situation  is  illustrated  in  Figure  I  where 
the  inset  gives  our  computed  B  values  for  years  when  the 
sunspot  number  was  nearly  constant  near  40  on  the  rising 
part  of  the  cycle,  away  from  recurrent  high-speed  streams 
and  with  no  lingering  flux  from  the  maximum  phase,  for 
comparison  with  the  main  part  of  the  figure  adapicd  from 
LSW99,  Figure  8  of  SCQ5  showed  that  71%  of  the  variation 
of  B  is  explained  by  the  varying  sunspot  number.  For 
constant  sunspot  number  wc  see  no  increase  in  the  coronal 
magnetic  field,  while  LSW99  do.  Nothing  in  their  comment 
on  our  paper  leads  us  to  doubt  our  technique  or  result  In 
fact,  it  has  given  us  the  opportunity  to  raise  a  serious 
question  about  their  methodology,  specifically,  their  use  of 
the  Sargent  recurrence  index  to  reconstruct  solar  wind  speed 
(combined  with  an  error  in  the  calibration  of  the  aa  index 
prior  to  1957  [Svalgaard  et  aL,  2004;  Jarvis,  2005; 
M  Lockwood  et  ai..  The  long-term  drift  in  geomagnetic 
activity:  calibration  of  the  aa  index  using  data  from  a 
variety  of  magnetometer  stations,  submitted  u>  Aimaks 
Geophysicae,  2006,  L.  Svalgaard  and  E,  W.  Clivcr,  Long¬ 
term  variation  of  geomagnetic  activity  (the  JHV-index)  and 
its  use  in  deriving  solar  wind  speed  since  1882,  manuscript  in 
preparation,  2006]),  which,  in  our  opinion,  is  the  likely  cause 
of  the  major  increase  in  the  infeircd  coronal  Odd  that  LSW99 
report,  but  which  we  cannot  confirm. 

[ia\  Both  our  method  and  that  of  LSW99  use  extrapola¬ 
tions  of  correlations  outside  the  range  for  which  they  were 
derived  and  thus  both  must  be  regarded  with  caution.  For 
this  reason,  in  SC05  we  substantiated  our  /Debased  results 
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I  igure  9.  Sargent’s  recurrence  index,  /,  calculated  as  the  13-rotation  running  mean  (heavy  black  line)  of 
Ihe  cross  correlation  for  each  pair  of  Bartels  rotations  since  1868  calculated  from  12-hourly  averages  of 
the  an  index  (grey  line).  The  value  for  2003  is  indicated  by  a  red  circle.  Note  that  the  '‘base  level”  of  /  is 
significantly  higher  before  ^1925  than  after  that  time.  This  has  the  effect  of  systematically  decreasing  the 
value  of  B  calculated  using  the  method  of  LSW99  before  1925  (compare  Figure  8), 


by  using  an  independent  method  based  on  magnetic  field 
measurements  in  the  terrestrial  polar  caps,  some  of  which 
were  made  during  the  crucial  interval  early  in  the  last 
century  for  which  LSW99  inferred  very  low  coronal  fields. 
Our  debate  with  Lockwood  and  colleagues  on  the  long-term 
evolution  of  the  coronal  magnetic  field  and  the  solar  wind 
may  he  resolved  within  the  next  few  years  if  our  prediction 
[Svalgaard  et  aLt  2005]  of  a  solar  maximum  with  peak 
sunspot  number  comparable  to  that  of  cycle  14  bears  out.  If 
so,  direct  measurements  of  solar  wind  properties  during 
conditions  similar  to  those  during  the  previous  minimum  of 
the  Glcissbcrg  cycle  would  take  the  estimates  of  IMF  B  out 
of  the  realm  of  extrapolation.  It  is  noteworthy  that  the  IDV 
index  (and  thus  B ,  regardless  of  regression  method)  for 
2006  (based  on  the  first  7  months  only,  but  expected  to  fall 
further  as  we  approach  solar  minimum)  ts  already  the  lowest 
in  the  last  94  yeans, 
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